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The geodesies of massless particles produced in collisions near a rotating black hole are solved 
numerically and a Monte Carlo integration of the momentum distribution of the massless particles is 
performed to calculate the fraction that escape the black hole to infinity. A distribution of in falling 
dark matter particles, which are assumed to annihilate to massless particles, is considered and an 
estimate of the emergent flux from the collisions is made. The energy spectrum of the emergent 
particles is found to contain two Lorentz shifted peaks centred on the mass of the dark matter. 
The separation of the peaks is found to depend on the density profile of the dark matter and could 
provide information about the size of the annihilation plateau around a black hole and the mass of 
the dark matter particle. 
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INTRODUCTION 



Intermediate mass black holes may be associated with 
a density spike of dark matter [Jj that can provide an en- 
hancement in the annihilation rate of dark matter parti- 
cles to gamma rays or neutrinos. The indirect detection 
of dark matter through its annihilations could provide 
information about the nature of the dark matter in the 
galaxy and its distribution. The physics of particle col- 
lisions in the gravitational field of a black hole has been 
extensively studied, see e.g. 0I3]IllIS]IS]I7][H][i] 
[To] . In the case of rotating black holes, collisions be- 
tween dark matter particles can yield large centre of mass 
energies [Hj and inside the ergosphere particles produced 
via the Penrose process could carry high energies fB] . In 
order to make use of any such collisions around a black 
hole it will be important to understand the fate of the 
particles produced, in particular the fraction of particles 
produced that escape the black hole compared to those 
which cross the event horizon. The escape fraction can 
be defined as the fraction of particles that escape having 
been produced in a collision of two dark matter parti- 
cles. The escape fraction depends on the momenta of the 
colliding dark matter particles and the distance from the 
black hole at which they collide. The escape fraction has 
been calculated for maximally rotating black holes where 
massless particles are emitted in the equatorial plane of 
the black hole [H] and is also known for Schwarzschild 
black holes [T^]. However a general analytical solution 
is not known for rotating black holes where the particles 
are not restricted to the equatorial plane. 

The goal of this paper is to find a numerical result for 
the escape fraction of massless particles produced in col- 
lisions around a rotating black hole and to investigate 
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the energy spectrum of any emergent flux. We assume 
an interaction of the form xx ~^ where x is a massless 
particle that could be a photon or a massless neutrino. 
(The analysis can also be used for neutrinos with small 
masses as the change in the result is negligible provided 
the mass of the dark matter is much larger than the neu- 
trino mass.) To simplify the analysis, we assume that 
the final state particles are emitted isotropically in the 
centre of mass frame of the collision. 

The paper is organised as follows: in Section |ll] the 
method of numerically solving the geodesic equation for 
a massless particle in a Kerr metric with arbitrary ro- 



tation is presented. Section [III] uses this to obtain the 
escape fraction for massless particles and introduces the 

effect of a boosted frame of reference due to the momen- 
tum of the colliding particles. In Section IV the dis- 



tribution of the momenta of the colliding dark matter 
particles is outlined. Section |V] presents the spectrum of 
emergent particles from the collisions and finds that the 
spectrum generally contains two Lorentz shifted peaks 
centred around the dark matter mass. In Section IVII we 
present our conclusions. 



II. NUMERICAL SOLUTION OF THE 
GEODESIC EQUATION 



To know whether a particular particle escapes the 
black hole or not, its geodesic equation must be solved 
and the behaviour as t — > oo found. This is most easily 
performed using Boyer-Lindquist coordinates as defined 
in the analysis of Ref. [13], which is followed closely here. 
In these coordinates the metric for a Kerr black hole has 
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where M is the mass of the black hole, a is the angular 
momentum per unit mass of the black hole, is the polar 
angle from the axis of rotation, (f) is the azimuthal angle 
and the following functions are defined |13) . 
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For simplicity the mass of the black hole is set to M = 
1, as a result a ranges between and 1 in units where 
^'^^"'•^ = 1 such that the Schwarzchild radius = 2. The 
trajectories of a particle are then defined by 3 constants of 
motion which are conserved, the forms of these quantities 
are [131 
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Here /i is the mass of the particle, E is the energy, L 
is the component of the angular momentum parallel to 
symmetry axis of the black hole and Q + is the total 
angular momentum squared when a = 0. Q characterises 
the motion in the 9 direction and for Q = particles in 
the equatorial plane will remain restricted to that plane 
|13j . The equations of motion for the particles can now 
be written in terms of A, which will become an affine 
parameter for massless particles, and are written here in 
the form given by Ref. |13) . 
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The final state particles of the collisions are assumed to 
be massless allowing the equations to be simplified by 
setting /X = (For massive particles the approximation 
/i = is good provided fj, <^ E and E is generally of the 
order of the dark matter mass. The results therefore re- 
main unchanged if a small mass is introduced.). In order 
to calculate the escape function numerically the equa- 
tions of motion for r and 9 were integrated numerically 
until the particle either crossed the horizon or escaped 
to some large value of r. The algorithm employed was 
based on an embedded adaptive step size Runge-Kutta 
formula as outlined in Ref. [14j. 



The integration was ended whenever the massless par- 
ticle crossed the horizon or achieved some maximum dis- 
tance from the black hole, for example r^^ax ~ 100 for 
collisions close to the horizon. This was increased for col- 
lisions further from the black hole. In the case of bound 
orbits the integration was ended after a suitable number 
of turning points in the r solution or a maximum value 
for A was reached. In practise the parameters of the algo- 
rithm, the minimum step size, A 



max ; ' max 



and the value 

of ^ that indicates a turning point were varied to find 
optimal values that did not introduce significant errors in 
the final result. This was checked against test cases and 
known solutions for the escape fraction ^2] [IT] shown in 
Figures. [T] and [2| 



III. THE ESCAPE FRACTION 

This method for solving the geodesic equation for the 
emergent particle then allows the escape fraction of mass- 
less particles produced in collisions around the black hole 
to be found. The escape fraction was calculated by a 
Monte Carlo integration over the final momentum of the 
produced massless particle pair. In the centre of mass 
frame of the colliding particles the massless particles are 
emitted back to back and carry energy equal to one half 
of the centre of mass collision energy. The momentum 
integration is therefore reduced to an integration over 
the direction of the massless particle, which we define 
in terms of two angles, 9 the angle from the radial axis 
parallel to the polar angle of the black hole coordinate 
system and (j) the azimuthal angle about the radial axis. 

For the integration only one of the pair of massless par- 
ticles was considered and a flat distribution of its angles 
in the centre of mass frame was generated. For each gen- 
erated direction the geodesic equation was solved numer- 
ically as above and the massless particle was recorded as 
either escaping or becoming trapped by the black hole. 
The escape fraction was then taken as the number of 
massless particles recorded as escaping divided by the 
total number generated. 

The uncertainty in the escape fraction due to the sta- 
tistical effects of the Monte Carlo integration was found 
by separating the massless particles into discrete bins, 
the escape fraction for each bin was determined then the 
average escape fraction taken. The variance in the escape 
fraction for each bin then gave an estimate for the uncer- 
tainty in the averaged escape fraction. This uncertainty 
does not account for any error in the numerical solution 
for the geodesic equations but since the escape fraction 
only depends on the general behaviour of each massless 
particle this was neglected. 

To test the implementation of the numerical method 
the escape function can be calculated for the a = 
Schwarzschild case. Analytically the conditions for es- 
cape are known, and can be converted to an escape frac- 
tion as a function of the radial position [T^]. Fig. [l] shows 
the behaviour of the escape fraction when the massless 
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FIG. 1: The escape fraction for a = 0, where N is the number 
of Monte Carlo iterations. The sohd line shows the escape 
fraction constructed from the analytical conditions of escape 



particle is emitted in an unboosted frame so that the 
emission is isotropic in the lab frame. The escape fraction 
matches the function determined from \T2\ with some 
variation due to the Monte Carlo integration. The ef- 
fect of the number of Monte Carlo iterations on the error 
in the escape fraction can be seen and gives an indica- 
tion of the number of iterations required to give a good 
agreement with the analytical solution. 



A. Locally non-rotating and boosted frames 

Having tested the machinery for the non-rotating case 
results can be obtained for black holes with a > 0. For 
a = 1 it will be useful to consider two cases for the nu- 
merical calculation, firstly where the massless particles 
are generated with no component of their motion outside 
the equatorial plane {Q = 0) and secondly where they 
are generated in all directions to compare with the result 
found in fTT'. For a rotating black hole there is another 
complication. So far the massless particles have been as- 
sumed to be emitted in a locally flat frame that is at rest 
in the Boyer-Lindquist coordinate system. For a rotat- 
ing black hole the corresponding frame is that of a locally 
non rotating frame (LNRF) |T3]. A LNRF represents a 
set of Minkowski coordinates such that the physics in this 
local frame can be described by special relativity. This is 
further complicated when the centre of mass frame of the 
collision is Lorentz boosted relative to the LNRF. In this 
case it is easiest to generate the massless particle momen- 
tum in the centre of mass frame where it is assumed to 
be isotropic and apply a Lorentz transformation to attain 
the momentum in the LNRF. 

The transformation from the isotropic distribution of 
momentum in the centre of mass frame to the distribution 
of the constants of motion in the Boyer-Lindquist coor- 
dinate system depends on the coordinates at which the 
collision takes place, rp and and A(Pi, P2) the Lorentz 
transformation between the centre of mass frame and the 
LNRF. Where Pi and P2 are defined as the momenta of 



the colliding particles in the LNRF. For details on trans- 
forming between a LNRF and a Boyer-Lindquist system 
see Ref. [13], while the method of calculating A(Pi,P2) 
we follow that in Ref. [Tl . 

The general procedure for the calculation of the mass- 
less particle momentum in the Boyer-Lindquist coordi- 
nates was carried out as follows. Firstly the momenta 
of the incoming particles were defined in the Boyer- 
Lindquist frame by setting the constants of motion L and 
Q then evaluating the equations for their momentum at 
the coordinates that the collision was assumed to take 
place. The incoming momenta were transformed to the 
LNRF as in Ref. [13, . With the two momenta of the in- 
coming particles defined in this locally flat space the rel- 
ative momentum of the centre of mass frame was found 
and the boost between these two frames determined. The 
massless particle direction was then randomly generated 
in the centre of mass frame, the boost between the LNRF 
and the centre of mass frame was inverted to allow the 
massless particle momentum in the LNRF to be found. 
Finally the massless particle momentum in the Boyer- 
Lindquist coordinates was calculated giving the constants 
of motion for the massless particle. The geodesic equa- 
tion for the massless particle could then be integrated 
numerically as before. 

The momentum of the incoming particles was deter- 
mined following the procedure in Ref. [llj which calcu- 
lated the velocities of the incoming particles that were 
restricted to the equatorial plane. The momenta of the 
two incoming particles was then used to produce the 4 
dimensional Lorentz boost between the LNRF and centre 
of mass frame. Fig.[2]shows the escape fraction as a func- 
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FIG. 2: Escape fraction for a=l with momenta of colliding 
particles fixed: Analytic solution in equatorial plane (Black 
line) [llj . Numerical result in equatorial plane (blue circles) 
and full numerical result allowing propagation in all directions 
(red squares). 

tion of r for collisions between two massive particles col- 
liding in the equatorial plane of a black hole with a — 1. 
The colliding particles have angular momentum Li = 2 
and L2 = —2 which correspond to the constant of mo- 
tion L and both have Q = 0. The escape fraction found 
through an analytical calculation [11] when restricted to 
the equatorial plane is also shown. The numerical result 
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matches the exact solution when the particles are con- 
fined to the equatorial plane. The numerical result when 
the massless particles are emitted in all directions gives 
a value for the escape fraction which is comparable in 
size to the result in equatorial plane, suggesting that the 
approximation used in Ref. dlj is reasonable. 





IV. DISTRIBUTION OF MOMENTA OF THE 2 4 6 8 10 12 14 16 

COLLIDING PARTICLES Q 



The previous section considered massless particles pro- 
duced in collisions of massive particles with fixed incom- 
ing momenta however the escape fraction depends on the 
Lorentz boost from the LNRF to the centre of mass frame 
given above, which is derived from the constants of mo- 
tion E, L and Q of the incoming particles. These con- 
stants are not fixed in general and the momenta of the 
colliding particles will have some distribution. This dis- 
tribution will also affect the centre of mass energy of the 
collision and the energy of the produced massless parti- 
cles as well as the direction in which they are emitted. 
In order to calculate the energy distribution of escaping 
particles and the general escape fraction it will therefore 
be important to consider incoming particles with some 
distribution of momenta. 

The form of this distribution is not known so as a first 
approximation we assume a flat distribution in the con- 
stants L and Q. The allowed range of L and Q was fixed 
by considering only particles with geodesies that would 
allow the particle to fall from some large r to the radius 
To at which the collision takes place. The constant of 
motion E was set to the mass of the dark matter which 
in turn is set to 1. The range of allowed constants of 
motion was calculated by fixing 9 to the equatorial plane 
then calculating ^ as a function of r starting at r ^ tq 
then running down in small steps of r to the value rg, 
if 1^ becomes very small then we assume that a turning 
point is reached and that a particle could not fall from 
infinity to tq. The distribution of allowed momenta was 
sampled by Monte Carlo integration. Figure. |3] shows 
the allowed range in L and Q for collisions occurring at 
rg = 1.1 as an example. This is characterised (for a = 1) 
by a maximum value of Q which occurs for L = — 1 as 
well as a maximum and minimum value of L with the 
largest range corresponding to Q = 0. Figure. |4] shows 
the escape fraction when the momenta of the colliding 
particles is averaged over all allowed values of L and Q. 
The escape fraction has only a weak dependence on the 
Lorentz boost and hence L and Q, of the colliding parti- 
cles so the full escape fraction is only a slight shift from 
the result found for a particular fixed boost shown in 
Figure. [2] 



FIG. 3: Allowed range in Q and L for incoming particles that 
reach ro — 1.1 with a = 1. The red line show Lmax and the 
green fine Lmin- 
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FIG. 4: Escape fraction for a — 1 with integration over mo- 
menta of the colliding particles including the effect of the 
Lorentz boost. 



V. SPECTRUM OF EMERGENT MASSLESS 
PARTICLES 

In order to calculate the energy spectrum of emergent 
particles we must integrate the effect of collisions over the 
distribution of dark matter surrounding the black hole. 
The flux arriving at some distance D, from the black hole 
reads [TT], 

"-A^^i P'(r)eir)dV, (6) 

where av is the cross-section for annihilation to mass- 
less particles which we assume to be energy independent, 
niy. is the mass of the dark matter particles, is the 
Schwarzschild radius, is the event horizon, e(r) is the 
escape fraction shown in Figure. |4] and p{r) is the dark 
matter density around the black hole. It is assumed that 
the escape fraction is spherically symmetric. The calcula- 
tion is performed for collisions occurring in the equatorial 
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plane only, to simplify the integration. By writing 
P{r) = PpiPair), 



(7) 



where poi^r) is a dimensionless function giving the shape 
of the density profile. The flux can be written in terms 
of a dimensionless integral containing all the of the de- 
pendence on the rotation of the black hole and the shape 
of the density profile as [llj . 



I(a,ri) 



r pQ{r) e(r) dr. 
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(9) 



In order to calculate the spectrum of the emergent mass- 
less particles we need to know <&(£') which depends on 
the convolution of the escape fraction and spectrum of 
massless particles produced in the collisions. To find this 
we define p(r, Ei,E2) as the fraction of massless particles 
that escape with energy between Ei and E2 per collision. 
This was found numerically by separating the escaped 
massless particles into bins of energy and dividing by 
the total number of collisions. The integral can then be 
rewritten as 

<^>{Ei,E2)xI{a,n,Ei,E2)^ / pl{r)p{r,Ei,E2)dr, 

(10) 

where summing over all the energy bins gives the total 
flux. 

The density distribution of dark matter around black 
holes has been extensively studied for non-rotating black 
holes see [T] [TJ] [IB]. Considering the case of an in- 
termediate mass black hole the density distribution can 
be described as follows [T]; close to the black hole there 
is an annihilation plateau with constant density ppi ~ 
m-^/ [av t) where t is the formation time of the black hole. 
This holds out to a radius rcut from which point the den- 
sity falls off with a power law p (xr~^ . The distribution 
here is assumed to be spherically symmetric (For a ro- 
tating black hole this may no longer be true but serves 
as a first approximation.). The flux will be dominated 
by contributions from the annihilation plateau and den- 
sity spike since the escape fraction is only small for radii 
much smaller than the annihilation plateau. It is now 
clear that if Ppi given in Eq. [7] is taken to be the density 
of the annihilation plateau, p{){r) will take the form [1], 
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if r < Tcuf 
if r > Truf 



(11) 



To calculate the value of I for a particular set of black 
hole parameters the value of Vcut needs to be specifled. 
In general this depends on the mass of the black hole and 
on the annihilation cross-section of the dark matter. For 
the super massive black hole at the centre of the galaxy 
this was estimated to be Vcut ~ 4 x \Q~^pc from Ref. [15j 



which is Tcut ~ 137rs in terms of the Schwarzschild ra- 
dius. 

The exact value is not important since the shape of 
the spectrum remains distinctive over a large range of 
r^ut and a measurement of the spectrum should in prin- 
ciple allow rcut to be determined. To show this, the spec- 
trum produced for a number of different values of rcut are 
shown in Figure. [5j The energy is expressed in units of 
the dark matter mass. 

Figure. 5(a)| shows the spectrum for a small value of 
fcut = 5 (in units where = 2), and hence is dominated 
by collisions close to rg = 5. Here the spectrum is char- 
acterised by a single asymmetric peak below the dark 
matter mass. The peak itself is broad and the mass of 
the dark matter is found to be roughly equidistant from 
the peak and the edge of the shoulder in the distribution. 
This structure arises due to several effects. Firstly there 
is a balance between the centre of mass energy and the 
gravitational red shift as suggested in Ref. [2j , this is ex- 
pected from energy conservation considerations that the 
energy available to the collision as measured far from the 
black hole will be 2m^. 

The result is that the spectrum is not shifted to higher 
energies by the increasing centre of mass energy close 
to the black hole horizon or red shifted by increasing 
gravitational potential. There is a shift of the peak to 
lower energies that arises from the Lorentz boost, since 
the centre of mass frame is in general boosted towards the 
horizon, massless particles emitted away from the black 
hole are red shifted to lower energies. Particles emitted 
towards the horizon are correspondingly shifted to higher 
energies. The effect of this is to split the spectrum of 
escaping particles into two peaks, a red shifted peak with 
energy below the dark matter mass and a blue shifted 
peak with energy above the dark matter mass. 

In Figure. |5(a)| only one such peak is visible this is 
because particles emitted towards the horizon are much 
less likely to escape the black hole. There is still a slight 
shoulder to the spectrum just above the dark matter mass 
which corresponds to the remnant of the blue shifted 
peak after the escape fraction has been taken into ac- 
count. 

As rcut in Figure. [5 (b)| increases to 34, both peaks be- 
come clearly visible and the spectrum becomes less broad. 
The height difference between the peaks also decreases 
as they get closer together. The peaks in the spectrum 
move towards the dark matter mass since the Lorentz 
boost to the centre of mass frame becomes smaller. The 
size of the Lorentz boost decreases as the average radial 
momenta of the colliding particles decreases. The escape 
fraction becomes less important as r increases and the 
massless particles initially emitted towards the black hole 
can now escape in many cases. It is still more likely that 
a massless particles emitted away from the horizon will 
escape than one emitted towards it so the blue shifted 
peak is generally smaller than the red shifted one. It 
also becomes clear that the spectrum is indeed roughly 
symmetric around the dark matter mass which could be 



6 




0.6 0.7 0.8 0.9 1 
Energy 

{h)rcut = 34 



900000 
800000 
700000 
_ 600000 
500000 

LU 

- 400000 
300000 
200000 
100000 




1.3 




0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 
Energy 

{c)rcut = 274 



the massless particles produced in the colhsions will now 
escape. This can be seen in the spectrum by noting that 
the red shifted and blue shifted peaks are now very nearly 
the same height. Comparing the spectrum with that in 
Figure. 5(b) it can be seen that the red shifted and blue 
shifted peaks continue to move towards the dark matter 
mass as the effect of the Lorentz boost diminishes. 

The spectrum could potentially be used to infer some 
of the quantities in the system. The dark matter mass 
is indicated by the minimum between the peaks and the 
value of r^ut can be estimated from the separation of the 
peaks as a fraction of the dark matter mass. 

The value of Vcut can be estimated as follows; the dark 
matter mass must be found first by locating the minimum 
of the spectrum, the energy at which the first peak occurs 
as a fraction of the mass of the dark matter then decreases 
as a function of r„if. By plotting the peak location as 
a function of Vcut from simulation the expected value of 
fcut can be found for a given peak separation. It should 
be noted that the distribution is not exactly symmetric 
and that the higher energy peak is located closer to the 
dark matter mass than the lower energy peak. However 
both could be used to find an estimate of Vcut and the ra- 
tio of the height of the peaks could also be utilised in this 
way. One problem with this analysis is that it assumes 
that the angular momentum of the black hole is known 
and varying a shifts the separation of the peaks. However 
this effect is much smaller than changing r^ut but would 
still introduce uncertainty in the estimation. The ability 
to resolve the peaks depends on the total flux available 
and the energy resolution around the dark matter mass 
of a particular measurement. In terms of the annihila- 
tion plateau the separation of the peaks was found up 
to r^ut — 10^ with a separation of A^; ^ O.OOlm^. For 
rcut — lO'^ a separation of ~ 0.03m^ was found. In 
the case that the peaks can not be resolved then the spec- 
trum can be approximated as a single peak with a width 
comparable to the peak separation centred at the dark 
matter mass. 

The total flux can be estimated by choosing some suit- 
able values for the parameters in Eq. [9j Setting the mass 
of the black hole to M = 40 x 10^ Mq, the annihilation 
cross-section of the dark matter as av = 10~^*cm^s~^, 
the distance from the black hole D — lOpc and the time- 
scale for the growth of the black hole as ~ 10^°years 



the flux can be written as [TTj, 



(12) 



FIG. 5: Spectrum of massless particles escaping to infinity 
with a cut off in dark matter density at various values of rcut ■ 
(a=l) 



useful in determining the mass from the spectrum. 

In Figure. |5(c) the spectrum is shown for rcut — 274 = 
137rs the collisions responsible for the majority of the 
flux are now occurring far from the horizon. The escape 
fraction becomes even less important and nearly all of 



where = 3.41km ^year ^. Integrating over r for 
^cut = 274 gives a total value of I w 2 x 10^ which 
gives a total flux of $ w 7 x lO'^km^^year^^ For lOTeV 
dark matter annihilating to high energy gamma rays the 
HESS experiment could be sensitive to such a large flux 
however the energy resolution would be smaller than the 
separation of the peaks [17]. For smaller values of rcut 
resolution of the peaks becomes easier but the total flux is 
reduced. At rcut = 34 the separation is potentially large 
enough to be observed but Z w Iff^ requiring greater 
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sensitivity or observation time. 




FIG. 6: I{rcut = 34, £1,^2) for a = 

The shape of the spectrum can also be considered for a 
non-rotating black hole and for a reasonable value of Vcut 
gives a flux very similar to the Kerr case. Figure. [6] shows 
the spectrum for a non-rotating black hole with r^Mt — 34 
and differs from Figure. 5(b) by having a slightly larger 
flux and the peaks closer to E = 1. This is due to the 
smaller Lorentz boosts associated with the Schwarzschild 
black hole and a faster growing escape fraction. This is 
found in contrast to the spectrum in Ref. [2] where the 
spectrum was found to have a narrow single peak. In 
considering the incoming dark matter particles to have 
momenta distributed over the whole range of allowed in 
falling geodesies the Lorentz shifting of the peak is found 
to have a splitting effect giving rise to two peaks that 
accounts for this difference. 



VI. CONCLUSIONS 

In conclusion we have shown that annihilation of dark 
matter to massless particles around a black hole can 
produce a distinctive signal independent of the particle 
model. The spectrum is found to be centred around the 
dark matter mass with a splitting into two peaks due to 
the Lorentz boosting of the centre of mass frame of the 
colliding particles. The shape of the spectrum retains its 
distinctive shape when integrated over a typical density 
profile for the dark matter and can in fact reveal details 
of the distribution, in particular the radius of the an- 
nihilation plateau thought to form around a black hole. 
The shape of the spectrum is not dependent on the total 
flux of massless particles from the annihilations which 
depends on several unknown parameters for a particu- 
lar source and is therefore a useful tool in discriminating 
such a signal from astrophysical backgrounds. 



In the analysis presented here the final state particles 
are assumed to be massless, however the final result is not 
changed for massive particles so long as the mass is small 
compared to the mass of the colliding particles. This 
will certainly be true for neutrinos in models where the 
dark matter mass is of the order of few GeV and above. 
The numerical calculation can be trivially expanded to 
include a mass term for the escaping particles and we 
find little change in the resulting spectrum even for final 
state particles with a mass half that of the dark matter 
particle. 

In calculating the energy spectrum we assumed that 
the annihilation cross-section of the dark matter was in- 
dependent of the centre of mass energy, for a particular 
dark matter model the energy dependence of the cross 
section could modify the shape of the spectrum. The 
cross-section would need to be found for each Monte- 
Carlo generated collision as the centre of mass energy 
changes for each pair of colliding particles. We leave this 
extension for future work. 

In the numerical analysis carried out the collisions were 
assumed to take place in the equatorial plane and the 
escape fraction and spectrum were taken to be indepen- 
dent of the initial polar angle 9q at which the collision 
takes place. The results were checked for collisions at 
60 = 0.6 for rcut — 34 and the spectrum retains the same 
structure with a lower total flux but the same separation 
between peaks. It would therefore appear that the ra- 
dial dependence of the energy spectrum dominates over 
any 9q dependence. For small r the spectrum will differ 
due to the dependence of the radius of the ergosphere 
on 6 however for reasonable value of r^ut the spectrum is 
dominated by collisions far beyond the ergosphere. 

The results also demonstrate that while large centre 
of mass energies are possible in collisions close to the 
horizon the emitted massless particles will in general be 
largely red shifted to energies around the dark matter 
mass. High energy massless particles were found in the 
spectrum with energies far larger than the dark matter 
mass due to the Penrose process, the non-zero escape 
fraction for these particles means that collisions around 
black holes could still be a source of energetic photons 
or neutrinos with energies several times the dark matter 
mass, but the flux of these will be small. 
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